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Abstract 

The  two-dimensional  continuous  wavelet  transform 
(CWT)  is  characterized  by  a  rotation  parameter,  in 
addition  to  the  usual  translations  and  dilations.  The 
CWT  has  been  interpreted  as  space-frequency  repre¬ 
sentation  of  two-dimensional  signals,  where  the  trans¬ 
lation  corresponds  to  the  position  variable,  and  the 
inverse  of  the  scale  and  the  rotation,  taken  together, 
correspond  to  the  spatial-frequency  variable.  The  in¬ 
tegral  of  the  CWT’s  squared  modulus,  with  respect 
to  all  variables,  gives  the  energy  of  the  original  sig¬ 
nal.  Therefore,  an  integration  on  a  subset  of  the  pa¬ 
rameters  gives  an  energy  density  in  the  remaining 
variables.  This  paper  deals  with  the  implementation 
of  the  two  basic  densities,  that  is,  the  position  (or 
aspect-angle)  and  scale-angle  densities. 

Introduction 

The  2D  continuous  wavelet  transform  (CWT)  has 
been  used  by  a  number  of  authors,  in  a  wide  vari¬ 
ety  of  physical  problems  [1].  It  has  also  been  ap- 
phed  in  designing  CWT-based  algorithms  for  detec¬ 
tion  and  recognition  of  targets  in  a  cluttered  envi¬ 
ronment  [2,  3,  4].  In  all  cases,  its  main  use  is  in 
the  analysis  of  images,  that  is,  the  detection  of  spe¬ 


cific  features.  Examples  of  features  are  hierarchical 
structures  or  particular  discontinuities  (e.g.,  edges, 
filaments,  contours,  and  boundaries  between  areas  of 
different  luminosity).  Of  course,  the  type  of  wavelet 
chosen  depends  on  the  precise  apphcation.  In  par¬ 
ticular,  the  detection  of  directions  requires  the  use 
of  an  oriented  wavelet  (Morlet),  whereas  an  isotropic 
wavelet  (Mexican  Hat)  suffices  for  pointwise  analysis. 

The  2D  CWT  is  a  representation  of  an  image  in 
a  feature  space  with  four  parameters:  scale,  orien¬ 
tation,  and  position.  Prom  the  energy-conservation 
theorem  [5,  6],  one  defines  various  energy  densities 
on  any  subset  of  the  four  variables. 

We  are  currently  investigating  these  energy  den¬ 
sities  in  order  to  built  a  CWT  2D  feature  detector 
for  a  given  target  type  in  various  conditions.  Those 
features  will  then  be  used  as  input  for  a  convolution 
neural  network  (CNN)  algorithm  [7].  These  densi¬ 
ties  will  also  be  used  to  design  an  ATR  algorithm  for 
detection,  classification,  and  recognition  of  targets  in 
FLIR  and  SAR  imagery. 

In  the  next  section,  we  recall  the  definition  of  the 
CWT,  and  give  its  energy  conservation  property.  We 
then  discuss  the  implementation  of  the  two  basic  en¬ 
ergy  densities,  and  the  complexity  of  the  related  al- 
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gorithm. 


Implementation  of  the  Two  Basic  Densities 


The  2D  Continuous  Wavelet  Transform 

The  2D  CWT  is  obtedned  by  projecting  an  image  s 
onto  a  family  of  dilated,  rotated,  and  translated  ver¬ 
sions  of  one  function  i/;,  called  the  analyzing  wavelet, 

5(a,0,6)  =  (1) 

=  — ^  (2) 

J  Cl 

(fk^^-^^*{ar-^ik))s{k)  ,(3) 

where: 


The  first  problem  one  faces  when  implementing  the 
CWT  is  that  of  visualization.  Indeed,  S{a,d,b)  is 
a  function  of  four  variables:  two  position  variables 
b  £  B?,  and  the  pair  {a, 9)  €  x  [0, 27r)  ~ 

One  may  say  that  the  CWT  has  rmfolded  the  signal 
from  two  to  four  dimensions;  this  feature  explains 
its  efficiency  in  decoupling  singularities,  but  at  the 
same  time  increases  the  memory  requirements.  As  a 
consequence,  some  of  the  variables  must  be  fixed  for 
visualization  of  the  CWT.  There  are  many  possibih- 
ties,  but  the  interpretation  of  the  parameter  space  as 
phase  space  given  above  suggests  two  natural  ways  of 
presenting  the  CWT,  using  two-dimensional  sections 
of  the  parameter  space: 


•  The  parameters  a,  9,  and  b  correspond  to  scale, 
rotation,  and  translation,  respectively.  The  hat 
stands  for  the  Fourier  Transform,  the  *  stands  for 
the  complex  conjugation,  and  r®  is  the  standard 
2x2  rotation  matrix. 


•  Ip  verifies  the  so-called  admissibihty  condition, 

=  {2-nf  f  <°°-  (4) 


In  practice,  we  will  choose  =  1. 

The  wavelet  transform  satisfies  the  following  en¬ 
ergy  conservation  property. 


=  [  (f]^\s(k)f  (5) 

=  jjj  ^d9(fb\S{a, 9,^)1“^  ■ 


This  leads  to  the  interpretation  of  |5(a, 0,6)|^  as  an 
energy  density  of  the  signal  s  in  position,  scale,  and 
orientation  variables.  There  is  a  one-to-one  map 
between  the  firequency  space  k  =  {kx,ky)  and  the 
scale-angle  space  {a,  9)  by  the  identification  |fe|  = 
a~^  and  tan~^(^)  =  9.  Therefore,  |5(a,6,6)p  = 

15(0“^,  0,6)1^  is  a  space-frequency  energy  density.  In 
other  words,  the  CWT  may  be  interpreted  as  a  phase 
space  representation  of  the  signal. 

It  is  clear  that  a  partial  integration  of  the  CWT 
energy  density  over  any  subset  of  the  variables  gives 
an  energy  density  in  the  remaining  variables.  Thus, 
there  are  four  ID  densities,  six  2D  densities,  and  four 
3D  densities. 

Among  these,  there  are  two  basic  densities:  the  po¬ 
sition  (or  range-aspect)  density  and-  scale-angle  den¬ 
sity.  In  the  following  section,  we  discuss  their  imple¬ 
mentation,  and  the  complexity  of  the  related  compu¬ 
tation. 


(i)  the  position  representation,  where  a  and  9  are 
fixed  and  the  CWT  is  considered  as  a  function 
of  position  b  alone.  The  corresponding  energy 
density  is 

E3iibx,by)=  /  -jd9\S{a,9,bx,by)\-^  ,  (6) 

Jo  Jo  ® 

(ii)  the  scale-angle  representation:  for  fixed  b,  the 
CWT  is  considered  as  a  function  of  scale  and 
angle  {a,  9),  i.e.  of  spatial  frequency 

Ei2{a,9)=  j  dH\S{a,9,bx,by)\^  .  (7) 

The  position  representation  is  the  standard  one, 
and  it  is  useful  for  the  general  purposes  of  image  pro¬ 
cessing:  detection  of  position,  shape  and  contours  of 
objects;  pattern  recognition;  image  filtering  by  resyn¬ 
thesis  after  elimination  of  unwanted  features  (for  in¬ 
stance,  noise).  The  scale-angle  representation  will  be 
particularly  interesting  whenever  scaling  behavior  (as 
in  fractals)  or  angular  selection  is  important,  in  par¬ 
ticular  when  directional  wavelets  are  used.  In  fact, 
both  representations  are  needed  for  a  full  understand¬ 
ing  of  the  properties  of  the  CWT  in  all  four  variables. 

Position  (or  Range- Aspect)  Energy  Density 

The  computation  of  E3i{bx,by)  requires  the  calcula¬ 
tion  of  S{aj,9j,bx,by)  for  j  =  l,....,iV,  that  is,  N 
two-dimensional  FFTs.  But  how  do  we  choose  the 
scale  aj  and  the  angle  9jl  To  answer  this  question 
one  needs  to  define  an  appropriate  sampling  in  the 
scale-angle  variables. 

Suppose  xp  is  a,  directional  wavelet,  that  is,  the  ef¬ 
fective  support  of  its  FT  xp  in  spatial  frequency  space 
is  contained  in  a  convex  cone  of  opening  angle  A(/3, 


with  apex  at  the  origin.  We  then  say  that  ‘tp  is  cen¬ 
tered  at  ko  if 

ko  =  J (fk  k\ip(k)\'^ .  (8) 

Its  width  in  the  x  and  y  directions  are  given  by  2wx, 
2wy,  where: 


(n  =  x,y)  . 

Taking  ko  =  (0,  ko)  along  the  y  axis  (see  figure  1),  we 
define  the  support  of  ^  as  the  ellipse 


Taking  intercepts  with  the  fc^-axis,  we  define  the  scale 
resolving  power  of  i/j  as 


For  a  Morlet  wavelet  of  anisotropy  e,  this  gives  [5,  6]: 

<“> 

ARP{iPm)  =  2mt~^(koVe)  .  (15) 

Using  the  angular-  and  scale-resolving  powers,  one 
obtains  a  tiling  of  the  frequency  space  by  the  wavelet 
V),  as  indicated  in  figure  2. 


Figure  2:  Tiling  of  frequency  space  using  the  parameters 
a  tind  9. 

Suppose  a  signal  s  ranges  in  spatial  frequency  be¬ 
tween  Pm  and  PM)  and  we  want  to  analyze  it  with 
a  wavelet  ip  with  scale  support  [pi,p2]-  The  extreme 
scales  that  are  needed  are,  respectively,  am  —  pm/ Pi 
and  om  =  Pm/P2)  and  the  scale  range  will  be: 

QM  _  Pm  / Pm  _  ^Psignal 

a^m  P2/  Pi  ^Pwavelet 

Therefore,  performing  the  discretization  on  a  dyadic 
scale,  the  signal  s  will  be  completely  analyzed  by  the 
family  of  wavelets  (‘filter  bank’)  {V’aj,e<(fc)})  where 

a,-  =  ^-2^  j  =  0,l,...,P-l,  (17) 

Pi 

dt  =  Ap-£,  ^  =  0,1,...,Q-1,  (18) 

with 


Figure  1:  Directional  wavelet  in  the  frequency  domain. 

As  for  the  angular  resolving  power  ARP  (ip)  =  Aip, 
we  consider  the  tangents  to  that  ellipse  and  obtain 

Jkl-wl 

ARP{ip)  =  2  cot-1  A -  (12) 

Wx 

~  2cot-^  — ,  for  ko^Wy.  (13) 

Wx 


P  =  integer  part  of  loga  ( 

\  ^Pwavelet  / 
27T 

Q  =  integer  part  of  — —  . 

Finally,  the  spatial  energy  density  is 


p-i  Q-i 


j=0  1=0 


(19) 

(20) 

(21) 


Scale- Angle  Density 

The  computation  of  Ei2{a,6)  requires  the  calcula¬ 
tion  of  S{a,9,bj)  for  j  —  As  the  geometry 

in  {hx,by)  plane  is  Cartesian,  the  sampling  will  be 
uniform,  that  is, 


jAa: ,  by,  =  jAy, 


(22) 


where  Aa;  and  Ay  are  the  horizontal  and  vertical 
bandwidths,  respectively,  of  the  wavelet  in  the  space 
domain.  We  are  still  faced  with  the  problem  of  calcu¬ 
lating  S(a,  6,  bj)  at  a  particular  point.  Let  us  consider 
the  wavelet  transform  at  o  =  bj 

S{aA  bj)  =  \  I  (f£r{r-efi))sr  (S)  ,  (23) 

CL  J  CL  ^ 

where  s-g,{S)  =  s(S  +  bj).  In  polar  coordinates,  we 
get 


n27r  n  n 

dpd(p  -  e)s^  {p,(p)  . 

a  a  ^ 


(24) 

(25) 

(26) 


Performing  the  change  of  variables 
p  =  e“  and  a  =  e^, 

we  obtain 
Sie\e,bj)  = 

/+00  |•21r 

/  du  dip  e“->*  (e“-^ ,  -  0)  e“  Sf,  (e“,  ip) . 

-ooJo  ’ 

Denoting 

G(u,  p)  =  p),  F{u,  p]  bj)  =  e“sj,  (e“,  <p),  (27) 

Eq.(26)  reduces  to 

Sie\e,bj)=  (28) 

/4-O0  r27r 

/  dudpG(u  —  v,p  —  6)F{u,p;bj)  . 

'OO  Jo 


In  this  way,  we  have  reduced  the  problem  of  comput¬ 
ing  S{a,6,bj)  into  that  of  computing  a  convolution. 
The  related  complexity  is  N  FFTs. 

Conclusion 

The  calculation  of  the  energy  densities  F12  (a,  9)  and 
F34{i)  uses  a  certain  number  N  of  2D  FFT.  This 
number  N  depends  on  the  support  of  the  wavelet 
and  signal  in  the  space  domain  for  Fi2ia,9),  and  the 
support  of  the  wavelet  and  the  signal  in  the  spatial 
frequency  domain  for  the  density  £34(6). 
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